Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Resolving the Galactic X-ray background 



M. Revnivtsev^'^, A. Vikhlinin^''^, and S. Sazonov 



1,2 



o 
o 

(N 
o 

CD 

Q 



(N 
> 
(N 

in 



o 



^ Mclx-Planck-Institute fiir Astrophysik, Karl-Schwarzschild-Str. 1, D-85740 Garching bei Miinchen, Germany, 
^ Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, 117997 Moscow, Russia 
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



Abstract. We use Chandra deep observations of the Galactic Center (GC) region to improve the constraints on 
the unresolved fraction of the Galactic X-ray background (also known as the Galactic ridge X-ray emission). We 
emphasize the importance of correcting the measured source counts at low fluxes for bias associated with Poisson 
noise. We find that at distances of 2'-4' from Sgr A* at least ~ 40% of the total X-ray emission in the energy band 
4-8 keV originates from point sources with luminosities I/2-10 kcv > 10^^ erg s~^. From a comparison of the source 
number-flux function in the GC region with the known luminosity function of faint X-ray sources in the Solar 
vicinity, we infer that Chandra has already resolved a large fraction of the cumulative contribution of cataclysmic 
variables to the total X-ray flux from the GC region. This comparison further indicates that most of the yet 
unresolved ~ 60% of the X-ray flux from the GC region is likely produced by weak cataclysmic variables and 



coronally active stars with L2-10 keV < 10 
is produced by discrete sources. 



erg s . We conclude that the bulk of the Galactic X-ray background 
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1. Introduction 

One of the largest extended features of the X-ray sky is 
the Galactic X-ray background, often referred to as the 
Galactic ridge X-ray e mission (GRXE hereafter), discov- 
ered in the late 1970's (ICooke. GrifBths. k PounS Il970l: 



Bleach et all 11973 iWorrall et al 



198^. The GRXE ex- 
tends over more than 100 degrees alo ng the Galactic plane 
but only a few degrees acro ss it (e.g., Warwick et al. 1985t 
Yamauchi fc Kovama 19931 ). 

The origin of the GRXE is the topic of a long-standing 
debate. The GRXE spectrum resembles the optically thin 
emission of thermal plasma with a temperature of 5- 
10 keV (e.g. lKovama et al.lll98dll989l ). However, such hot 
plasma cannot be gravitationally bound to the Galaxy, as 
suggested by strong concentration of the GRXE towards 
the Galactic disk and bulge. A number of models, still as- 
suming a diffuse o rigin of the GR XE, were proposed (see 
e.g. the review bv lTanakal liooi). but none of them was 
successful in explaining all of observed properties of the 
GRXE. 

An alternative explanation of the GRXE, proposed 
soon after its discovery, is that it is t he superposi- 
tion o f weak Galactic X-r a y sou rces 
1982t 



Worrall fc Marshalll 1 19831: iKovama et al 



Worrall et al 
19861: 



19921: iMukai fc Shiokawal Il993l) 



Ottmann fc Schmitt 
However, until recently there remained a large uncer- 



tainty with regard to the expected contributions of differ- 
ent classes of faint Galactic X-ray sources to the GRXE. 

( 2006I) . 



In recent work of . 
Revnivtsev. Molkov. fc Sazonovl 



Krivonos et al 



Revnivtsev et al 
( 200^ 



and 

based on X-ray observations 
INTEGRAL satellites, it was 



liooj^ 

the RXTE and 

that the GRXE closely traces the near- 
and consequently the stellar mass 
Furthermore, the observed 



from 

demonstrated 
infrared emission 
distribution in the Galaxy. 
GRXE to stellar mass ratio is compatible with the X-ray 
luminosity function of cataclysmic variables and coronally 
active stars in the vicinity of the Sun, measured with 
RXTE and ROSAT (jSazonov et al.ll2006l ). These findings 
therefore suggest that the bulk of the GRXE could 
i ndeed be produc ed by emission from descrete sources 



( Revnivtsev et al.ll200' 



A straightforward way to test the origin of the GRXE 
is to check which fraction can be resolved into discrete 
sources in deep Chandra images. The first such study 



was by lEbisawa et al. ( 2001 . 2005 ) who used a ~ 200 ks 



Send offprint requests to: mikej@mpa-garching.mpg.de 



Chandra observation of a Galactic plane region. They in- 
ferred that less than 15% of the GRXE in that region 
could be resolved into point sources with fluxes higher 
than (3-5)xl0~^^ erg s~^ cm~^. 

Perhaps the best place for studying the Galactic X-ray 
background is the Galactic Center (GC) region, because: 
1) the GRXE spectrum in this region is fairly typical (e.g., 
Tanakal[200^ . 2) the GRXE intensity is so high that the 
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Fig. 1. Chandra image of the Galactic Center region in 
the energy band 4-8 keV (efhciency corrected). The area 
of our study is outlined by the solid line. 

contribution of the extragalactic X-ray and the instrumen- 
tal backgrounds is very low, and 3) the high concentration 
of X-ray sources near the GC implies that most of them 
are within the Galactic nuclear stellar cluster and thus at 
the known distance from the Sun; this allows one to easily 
convert the source log N — log S relation to the intrinsic 
luminosity function. A disadvantage of the GC region for 
the GRXE studies is the presence of a large number of 
supernova r emnants and different types of no n-thermal 
phenomena (jMuno et al.ll2004 IPark et al.ll2004[ ). 

Chandra has accumulated a long exposure time in the 
several arcmin^ near the GC. Previo us analyses of the 



combined ~ 600 ksec worth of data ( Muno et al. l2003l 



20041 iPark et al.ll2004f) have shown that ~ 20-30% of the 
total X-ray flux can be resolved into point sources in a sub- 
region that is ma ximally free from supernova remnants 
(region "Close" in lMuno et al.l 120041 ). Since these works, 
an additional 330 ksec of GC observations have become 
available, which allows one to go deeper in resolving the 
Galactic X-ray background. In this paper, we analyze all 
publicly available Chandra data in the GC region in an 
attempt to resolve as much of the GXRE as possible. 

2. Data reduction 

We used Chandra ACIS-I observations of the GC re- 
gion with the aim point close to Sgr A* (observation 
ID #945, 1561, 2282, 2284, 2287, 2291, 2293, 2943, 
2951, 2952, 2953, 2954, 3392, 3393, 3663, 3665, 4500, 
4683, 4684, 5360, and 6113). The Chandra data were re- 
duced following a sta ndard procedure fully described in 
Vikhlinin et"aIT(|2005^ . The only difference is that the de- 
tector background was modeled using the stowed dataset 
(pittp:/ /cxc. harvard.edu/contrib/maxim/stowed). The to- 
tal clean exposure time is 918 ks. Figure [T] shows the com- 
bined image in the 4-8 keV energy band 4-8 keV. 




Fig. 2. Part of the image of the Galactic Center region in 
the energy band 4-8 keV, with the detected sources shown 
by circles (note that the source localizations are actually 
much more accurate than ~ 2.5 arcsec, the circle radius). 




Fig. 3. Comparison of source detections by wavelet de- 
composition (circles) and the CIAO task wavdetect 
(boxes) . Most of the sources are found by both detection 
algorithms. 



The GC re gion is known to be rich in superno va rem- 
nants (see e.g. lPark et al.|[20M : iMuno et al.ll2004l ). In or- 



der to minimize the contribution of X-ray emission from 
plasmas heated by the supernova remnants, in our anal- 
ysis we: 1) considered only a sub-region to south-west of 
the GC (see Fig. [1]), subtending 9.77 arcmin^, where the 
contribution of low-temperature plasmas is known to be 
sma ll from the weak ness of low-energy X-ray lines (see 
e.g. iPark et al.l |2004[ ) and 2) used only the energy band 
4-8 keV. 
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Fig. 4. Simulated image of the Galactic Center region (left) vs. the real image taken by Chandra (right). One can see 
a close similarity between them. The Galactic Center is in the upper left corner. 



To maximize the sensitivity to point sources, we used 
only the data within R < A' from the optical axis, where 
the Chandra angular resolution is within 1" (FWHM). 
Using a stacked image of bright (> 100 cnts) sources, we 
verified that the effective Point Spread Function (PSF) in 
the combined image can be well modeled as a Gaussian 
with a - 0.51". 

3. Source counting strategy 

To determine the total point source flux we cannot simply 
coadd the observed counts from detected sources, because 
of strong statistical biases discussed below. A better ap- 
proach is to integrate the reconstructed log N — log S func- 
tion down to the sensitivity limit. 

Because of its excellent angular resolution and low in- 
strumental background, Chandra can detect point sources 
yielding just a few total counts. Such low detection thresh- 
olds lead to strong biases in the derived log N — log S re- 
lated to the Poisson counting statistics. Generally, the raw 
number-flux relation near the detection threshold is sig- 
nificanly below the true log N — log S because of these ef- 
fects. This observational bias and methods for correctin; 



it have been extensively studied fe.g. lHasinger et al.lll99 
Vikhlinin et a l. 1995, and references therei n), and in par 



ticular for the Chandra image a nalysis in iMoretti et al.. 
(|2Q02[ ): lKenter fc Murravl (|2003h : Ifiauer et all (|2004[ l. An 



additional bias in our case arises from source confusion, 
non-neglibile in the GC region even with the Chandra^s 
angular resolution (Figl2]). 

Our procedure for recovering the true log N ~ log S 
function in the GC region is based on Monte-Carlo simula- 
tions, which allows us to derive accurate corrections for the 
statistical and source confusion biases. We start with sim- 
ulating the images containing realistic parent point source 
populations. We assume that the shape of the log A^— log 5 



relation is independent of distance from Sgr A* (which cor- 
responds to distance-independent luminosity function in 
this region) . The normalization of the parent log A — log S 
is assumed to vary as dN/dQ, oc R~^, consisted with 
the previous st udies of both the descret e source popu- 
Muno et a l!" '2003', ' 2006a[ l and the total 
Neronov et ah ,2005h in the GC region. 



lations (see e.g 
X-ray flux (e.g. 
The shape of the parent log A^ — log S is assumed to be a 
power law, dN/dS oc S~°' with a — 3.0, at bright fluxes 
(/ > 40 cntjil- At lower fluxes, we allowed for a nearly 
arbitrary shape of the intrinsic log A^ — log 5*. Namely, we 
assumed that S'^dN/dS — const in the flux intervals 40- 
10, 10~3, 3-1, and 1-0.3 cnts but the normalization within 
each interval is arbitrary. The normalizations of the par- 
ent log N — log S in each flux range were varied randomly 
in each realization. For each simulated source, the number 
of detected photons was drawn from the Poisson distribu- 
tion. If the total number of simulated photons was below 
the total observed intensity, the "missing" flux was added 
as a diffuse component with the surface brightness oc R^^. 
An example of the simulated image is shown in FigH) 

Our source detection is ba sed on the wavelet deco mpo- 
sition algorithm described in lVikhlinin et al~ ( 1998h . The 
detection threshold (specified as a required statistical sig- 
nificance) was chosen so that less than 1 spurious detection 
is allowed over the region of interest. Sources detected in 
the Chandra image of the GC region are shown by cir- 
cles in FiglSl We also checked that wavelet decomposition 
provides equal or better sensitivity to the point sources 
compared with the CIAO task wavdetect (Fig. [3]). 



^ Assuming the power law spectrum wit h the photon inde x 
r = 2 absorbed with iVn ~ 5 x 10^^ cm'^ (|Muno et al.ll2004l '). 
40 counts in the 4-8 keV band corresponds to an unabsorved 
flux of 7 X 10"^'^ erg s"^ cm"^ in the 2-10 keV band. 
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Fig. 5. Differential number-flux function of sources in the 
region of our study (points with error bars) along with 
a well-matching simulated number-flux distribution (solid 
line) and the corresponding parent (before the detection 
procedure) distribution of sources (dashed line) . Gray area 
denotes the manifold of trial shapes of differential lumi- 
nosity functions of simulated sources which we considered 
as satisfactorily describing the observed dN / dS function 



The fluxes of detected sources were measured within 
the 2.5" aperture (which should contain nearly 100% of 
the total source flux for the observed PSF width, see §[5]). 
The local background, provided by the largest scale of the 
wavelet decomposition, ^ 16", was subtracted from the 
source flux. 

The raw differential number-flux relation derived from 
the real observation is shown by crosses in FigO From our 
simulations, we can determine the "response" to which 
shape and normalization of the parent log N — log S best 
describes the observed flux distribution. These results are 
discussed below. 



4. Results 

The region outlined by the solid line in Fig. [T] contains 
in total 63.5 kcnts in the energy band 4-8 keV. The esti- 
mated particle background image of the same region con- 
tains 10.6 kcnts. Thus the net background-corrected num- 
ber of counts is 52.9 kcnts. The total flux from all the 
sources detected in the region is ~ 14.3 kcnts. 

4.1. Resolved fraction 

Figure [5] shows an example of a simulated number-flux 
function that is consistent with the distribution measured 
by Chandra. Also shown is the corresponding intrinsic 
log TV — logs' relation. By sampling over many such well- 
matching (the reduced difference from the observed 
number-flux function is less than 1.5) trial distributions, 
we determined the allowed range of the cumulative flux 



of "parent" sources in the GC region (shaded region in 
Fig. El). 

The outcome of this analysis is presented in Fig. [SI 
We conclude that at least 40% and possibly 100% of the 
total X-ray flux from the GC region is produced by point 
sources. We note that in obtaining this result we limited 
ourselves to sources with average intensities of more than 
3 integrated counts, because the number densities of yet 
weaker sources are poorly constrained by the available 
data. Therefore 40% is actually a conservative lower limit 
on the contribution of point sources to the total X-ray flux. 
Below we discuss the expected contribution of sources with 
luminosities below the Chandra detection threshold. 



4.2. Luminosity function of point X-ray sources 

We now address the luminosity function of faint X-ray 
sources in the GC region. Since the volume density of 
stars within several tens of parsecs of Sgr A* is several 
orders of magnitude higher than in other places along the 
line of sight in that direction, we can safely assume that 
the absolute majority of stars and X-ray sources in the 
sky region of our study (located at a ~ 10 pc projected 
distance from Sgr A*, see Fig. [TJ physically belong to ei- 
ther the nuclear stellar cluster (NSC) or nuclear stellar 
disk ( NSD) components of the Galaxv (iGenzel &: Towned 



19871 : iLaunhardt. Zvlka. k Mezged 12002 )7 This allows us 
to readily estimate the stellar mass contained in the vol- 
ume of the Galaxy covered by our observations and con- 
sequently the luminosity function of X-ray sources in the 
GC region normalized by stellar mass. 

For the NSC we adopted the following density profile: 

Pr 

Pmc 



1 + {r/rcT 

where pc = 3.3 x lO^A/g/pc^ and = 0.22 pc. The slope 
was assumed to be n = 2 at r < 6 pc and n = 3 at larger 
distances. The total mass of the NSC within 200 pc of Sgr 
A* is thus 6 X IO^Mq. T his value is actually uncertain by 
a factor of ~ 2 ( see e.g. Lindavist. Habing. fc Winnberd 



19921: ILaunhardt. Zvlka. fc Mezgeill2002l ) 



The NSD was assumed to have the density distribution 
(r and z are measured in parsecs) 



Pnsd = Pdr 



-\z\/zd 



where pd — 3OOM0/pc'^ and Zd — 45 pc. At r < 120 pc, 
the slope a = 0.1, at 120 pc < r < 220 pc, a — 3.5, 
and at r > 220 pc, a = 10. The total mass of the NSD 
is thus 1.4 X 10^M f7|. In reality this quantity is uncer- 
tain by some 50% (ILaunhardt. Zvlka. &: Mezgeil |200^ . 
We took the distance to the Galactic Center to be 7.6 



kpc (jEisenhauer et al.l 120051 ). 

The adopted mass model gives a total mass of stars 
enclosed in the region of our study of ~ 4.4 x IO^Mq. 
This value is uncertain by a factor of ~ 2 mainly due to the 
uncertainty in the NSC mass. We note that a similar mass 
model was found to provide a satisfactory description of 



Revnivtsev et al.: Resolving the Galactic X-ray background 



5 



X 

«^ 
"(5 

o 

+-» 
>-(— 
o 
c 
g 
*+-» 
u 




10 100 

Sources flux limit, counts 

Fig. 6. Cumulative flux (relative to the total flux) from 
the detected sources in the real and simulated observa- 
tions and the allowed range for the parent source distri- 
bution. Note that one count during the whole observa- 
tion corresponds to an intrinsic (corrected for the line-of- 
sight extinction of A'^h = 5 x 10^^ cm^^) source luminosity 
~ 1.3 X 10^° erg s^^. 



the surface densit y distribution of de tected X-ray sources 



in the GC region (jMuno et al.ll2006al ) 



Using this estimate of the enclosed stellar mass and the 
Chandra source number-flux function, we obtained from 
simulations the allowed range for the luminosity func- 
tion of X-ray sources in the GC region at luminosities 
-^^2-10 kov > 10^° erg s~^, which is shown in Fig. [71 It is 
important to note that the overall normalization of this 
luminosity function is uncertain by a factor of ^ 2 due to 
the uncertainty in the stellar mass enclosed in the studied 
region of the Galaxy. For comparison we show in Fig. [7] 
(points with error bars) the result of simply dividing the 
differential log A^-log S function measured by Chandra by 
the enclosed stellar mass. This naive determination is af- 
fected by the bias discussed in the previous section and 
thus underestimates the true luminosity function. On the 
other hand, the allowed range for the luminosity function 
inferred from the Chandra data for the dense GC region 
is remarkably compatible with the luminos ity function of 



faint X-ray sources in the Solar vicinity (jSazonov et al 



20061 Fig. [7]), taking into account the uncertainty in the 
NSC mass. 

The total absorption-corrected (assuming a line-of- 
sight absorption of A^h = 5 x 10^^ cm~^) X-ray flux in 
the energy band 2-10 keV from the studied GC region is 
Fx = (5.8 ± 0.4) X 10-12 erg s'^ cm'^. For the GC dis- 



tance of 7.6 kpc (Eiscnhaue r et al1l2005[ ). this corresponds 
to a total X-ray luminosity L2-10 keV = (4.0 ± 0.2) x 10'^^ 
erg s"i. Therefore the total X-ray emissivity in the stud- 
ied volume of the Galaxy is L^/M = (9.1 ± 4.6) x 10^'^ 
erg s~^Mq^, where we included the uncertainty in the 
NSC mass. This derived value of the X-ray emissivity 
per unit stellar mass agrees within the uncertainties with 
the cumulative emissivity of faint X-ray sources (cata- 
clysmic variables and coronally active stars) near the Sun: 
(4.5 ± 0.9) X 10^7 erg s-^Mg^ (ISazonov et al.ll2006^ . 

We should note here that in the studied GC region 
there may be a non-negligible contribution from warm dif- 
fuse plasma heated by supernova remnants. According to 
iMuno et all (f2004) , - 15% of the total 2-10 keV flux is 
probably due to a ~ 0.8 keV plasma. This would mean 
that the X-ray emissivity of point sources in the GC region 
is actually somewhat smaller, L^/M = (7.7 ± 3.9) x 10^'^ 
i. e. ev en closer to the value measured by 



erg s 



Sazonov et al. I (|2006h in the Solar neighborhood 
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Fig. 7. Luminosity function of weak X-ray sources in the 
Solar neighborhood (red area) in comparison with the al- 
lowed range for the luminosity function in the GC region 
(grey area). The points with error bars show a luminos- 
ity function estimated directly from the GC number-flux 
function, without correcting for the bias associated with 
Poisson counting noise at low fluxes. The normalization of 
the GC luminosity function per unit stellar mass is uncer- 
tain by a factor of ~ 2 due to the uncertainty in the mass 
of the nuclear stellar cluster. 



5. Discussion 

We showed that at least 40% of the total X-ray emis- 
sion from the Galactic Center region is produced by 
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point sources with luminosities L2-10 kcV > 10^^ erg s" 
The inferred lumin osity function of such sources (see also 
Muno et al.ll2006ah is compatible with t hat measured in 



the Solar vicinity ( Sazonov et al. 20061 ). Moreover, the 



data are consistent with the hypothesis that sources with 
luminosities below ^ 10^^ erg s^^, the effective Chandra 
detection limit, provide the rest of the total X-ray flux. 
Such sources are in fact expected to be present in the re- 
quired numbers in the GC region. 

According to studies in the Solar neighborhood 

( Sazonov et al. 20061 ). most of the X-ray sources with lu- 



minosities L2-10 keV < 10 erg s~ are coronally active 
stars. Late-type stars with convective envelopes can have 
coronae where th e plasma ca. n be heated up to X-ray 
temperatures (see Giidell 20041 for a review). Such stars 
are usually fast rotators, as is required for high coronal 
activity, because they are either members of binary sys- 



tems or relatively y oung (see e.g. IWalter &: Bowveijll981 
Queloz et al.lll998l ). The relative fraction of binary stars 



with active coronae near the Sun is quite high: at least 
every 80'th star has an X-ray (2-10 keV) luminosity in 
the range 10^"^ erg s~^ < ^2-10 koV < 2 x 10^° erg s~^. 
The relative fraction of single stars with active coronae is 
even larger. 

The relative contribution of coronally active stars to 
the total X-ray emission of a given volume of th e Galaxy 



is estimated to be 30-60% ( Sazonov et al. 20061 ). For this 



fraction to be smaller, the fraction of stars with convec- 
tive envelopes (usually late-type, low-mass stars) should 
somehow be reduced, whereas the population of low-mass 
stars at distanc es 4.5-8.7 pc fr om Sgr A* appears quite 
normal (e.g. .Philipp et al. 1999f) . 

As mentioned above, the luminosity function of GC 
sources agrees within the uncertainties with that measured 
in the Solar vicinity (Sazonov et al...200&) . We point out 
that because the vast majority of sources detected in the 
GC region are located at the same and known distance 
from us, a luminosity function determined there can be 
more accurate than a luminosity function determined else- 
where. 
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